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SUMMARY 

The entire plasma-sheath region is treated by a uniform method for an infi- 
nite plane electrode that adjoins a semi-infinite plasma in an electrical field 
normal to the electrode and a magnetic field parallel to the electrode. The ef- 
fect of collisions is neglected in calculating the velocity distribution func- 
tion. Both ions and electrons drift perpendicularly to the electric and magnetic 
fields, but there is no current to the electrode. A self-consistent potential is 
calculated, as well as the macroscopic properties of the electron and ion fluids, 
density, velocity, and pressure tensor components, for a hydrogen plasma in ther- 
mal equilibrium. An equilibrium velocity distribution and an isotropic pressure 
tensor are shown to be established at a distance jnore than several times the 
Larmor radius for each species; this minimal distance is reduced for ions in a 
sheath of large electron surplus. If the plasma density is sufficiently low, the 
entire transition region is electrically charged, whereas, near an anode in a 
plasma of higher density, the charge extends only to several times the Debye dis- 
tance. In the latter case, there is also an exterior region of charge neutrality 
and potential disturbance that extends to several times the electron Larmor ra- 
dius for a large anode potential and to several times the ion Larmor radius for 
other electrode potentials. 


INTRODUCTION 

When an electrode is placed in a neutral plasma, it tends to disturb the 
plasma in a way that depends on the electrode potential and the conditions in the 
plasma. Tonks and Langmuir (ref. l) showed that the disturbance region may be 
considered to be divided into a nonneutral region designated as the sheath and a 
neutral plasma region with a disturbance potential. The purpose of this report 
is to describe the structure of the steady-state condition of the entire distur- 
bance region adjacent to an infinite plane electrode when a magnetic field exists 
parallel to the electrode and when the density is sufficiently low that the ef- 
fect of collisions can be neglected. It is expected that some insight will be 
provided for the solution of the technically significant case where collisions 
and time variations occur. 

A number of studies have dealt with the electrostatic probe and with the 
sheath that determines the probe characteristics (see ref. 2 for an extensive 
bibliography and discussion). The anode probe in a magnetic field is treated by 



Bohm, Burhop, and Massey (ref. 3). This treatment is concerned only with the 
neutral region exterior to the sheath and is appropriate only when the ion Larmor 
radius is larger than the Debye radius. Bertotti (refs. 4 and 5) discusses the 
probe in a magnetic field without detailed examination of the sheath other than 
some general remarks , which are compared herein with the results of the present 
calculations. The transition region considered herein includes both the region 
of net electrical charge and the region of disturbance of electrical field by the 
electrode. 

The method used in this study is similar to that of Bernstein and Rabinowitz 
(ref. 6) in that the entire disturbance region generated by the electrode without 
a magnetic field is treated as a whole by means of the velocity distribution 
function and by ignoring the effect of collisions. Hall (ref. 2) points out that 
this method is suitable only for a probe of limited size, such as the spherical 
or cylindrical ones considered by Bernstein and Rabinowitz, but is unsuitable for 
a plane electrode, because a plane electrode that neutralizes impinging particles 
does not permit the return of attracted particles to the plasma, so that near the 
electrode the velocity distribution of the attracted species is one sided (no 
particles proceed away from the electrode). This situation will also persist in- 
definitely into a collisionless plasma since no mechanism is provided to effect 
the transition from the one-sided to the full velocity distribution. In the case 
of the spherical probe, the number of attracted particles that would strike the 
probe so decreases with distance that a transition to a full distribution is ac- 
complished without collisions. 

In case there is a magnetic force field directed parallel to a plane elec- 
trode (see fig. l), a transition to an equilibrium distribution in the plasma is 
possible without collisions because the motion of the particles normal to the 
electrode is inhibited. In this situation, however, no electrode current is 
possible; steady-state currents to the electrode are present only when there are 
collisions. The effect of a very small collision rate cannot be ignored alto- 
gether, because the effect of collisions is cumulative. If the collisions were 
ignored entirely, the plasma would remain frozen in whatever velocity distribu- 
tion might be initially imposed. If collisions occur so infrequently that the 
orbit of the average particle is slowly altered, the motion might be calculated 
over short periods of time and the slow changes ignored. After a long period of 
time, the velocity and spatial distributions will have altered substantially 
from the initial condition to an equilibrium distribution (except for the dis- 
turbance created by the electrode). 

The assumptions of the problem to be treated are the following: 

(1) The plasma and the sheath are at such a low density that collisions have 
a negligible effect on the velocity distribution at any instant of time. The 
long-period effect of collisions is accounted for by assuming an equilibrium dis- 
tribution far from the electrode. 

(2) Changes with time are sufficiently slow that at any instant of time 
there is approximately a steady-state condition. 

(3) An infinite plane electrode exists in the z,x plane, and the problem 
Is one dimensional. 
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(4) The electrode has absorbed or neutralized all particles that may have 
impinged upon it. 

(5) The plasma is neutral far from the electrode. 

(6) The plasma as a whole does not drift parallel to the magnetic field. 

(7) Thermal and drift velocities are small compared with the speed of light. 


DERIVATION OF DISTRIBUTION FUNCTION 


Equilibrium Distribution 

The electric field E is directed normally to the electrode in the plane 
y = 0, so that the electric potential cp is a function of y alone. The mag- 
netic field B is taken to be in the direction i z so that the magnetic poten- 
tial A satisfies 

I = A(y)i x 
B = B(y)i z 


dA 

ay 


i 


z 


(All symbols are defined in appendix A. ) Then, in Gaussian units , the Maxwell 
equations are 

— % = 4itec 2 (n e - (l) 

dy^ 



4jte(n e v x ^ e 


z i n i v x,i> 


(2) 


where 

cp(0) = A(0) = 0 
iHQ ^ = -Bq 

y^oo * 

lim — — Eq 

jr« y 

and e is the proton charge, c is the speed of light, n is the number density, 
Z is the number of charges per ion, v x is the average velocity in the x- 
direction, i x , i , and ± z are unit vectors, and the subscripts e and i 
refer to electrons and ions, respectively. When a particle of mass m moves in 
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a steady electric field, the energy 


£ = 



+ Zecp 


is constant. Of the generalized momenta, 


p x = mv x + ZeA 

p y = mv y > 

P z = mv z 


(3a) 


(3b) 


p x and p z are constants of motion. Then, for a given field, the velocity Vy 
is obtained from equation (3a) and the integrals of motion: 


= — (£• - Zecp) - 


(p x - ZeA) 2 + 


m 


m 


(4) 


For the steady, collisionless plasma the velocity distribution function f 
depends on the variables y, v x , Vy, and v z . (See appendix B for conditions 
under which collisions can be ignored. ) When these variables are transformed to 
the variables £, p x , Vy, and p z , the Boltzmann equation may be written 

Bf dS df df df dv y 

cig dt + dp x dt + dp z dt + dvy dt ~ 

Because £, p x , and p z are constants of motion, the solution is 

f = ) 


This solution indicates that the particles will remain permanently in any initial 
distribution. Actually, the collisions will cause a slow change to an equilib- 
rium distribution; therefore, out of all the possible distributions, it is only 
this one that is of interest here. For a subsystem of a plasma system, the sta- 
tistical distribution function is an exponential function of a linear combination 
of the additive integrals of motion (see ref. 7). For the present case, each 
particle may be considered to be a subsystem because of the assumed weak interac- 
tion between the particles; also, f is the statistical distribution function for 
a single particle, and the integrals of motion are £, p x , and p z . The velocity 
distribution function for particles in equilibrium and with a negligible colli- 
sion rate is therefore 

( e - ap x - pp z \ 

f = C exp(- — J 

where C, a, (3, and 0 are constants. When the constants are evaluated, 6 is 
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found to be kT, where k is the Boltzmann constant and T is the temperature, 
and a and 0 are the plasma drift velocities along the x and the z axes, 
respectively. 


Effect of Electrode 

The effect of the electrode is assumed to consist in the removal of all par- 
ticles that would collide with it. Thus, for the sheath, 


f = CH(p ,p ,£)exp 


a-P x + 3P 7 . - £ 


f 


e 


where H = 0 for particles that collide with the electrode and H = 1 other- 
wise. 


For any particle with assigned values of g, p x , and p z , the value of v 2 
may be found at any point by equation (4); if v 2 < 0, the particle cannot reach 
the point. Admissible particles are limited to those that cannot reach the elec- 
trode, that is, those for which 


0 5 v y(o) - I e - ± {4 + 

The electrode effect is then obtained by identifying 
function: 


p D 

H(U) 


with the Heaviside 


where 


H(U) = 0 

for 

U < 

H(U) = 1 

for 

U > 

2 , 2 



P x + P z 

U s x z 

- 2 

g 

m 2 


m 



(5) 


The situation is actually somewhat more complex than condition (5) would in- 
dicate, and consequently the present formulation has a restricted range of appli- 
cation. The limitation on the class of fields for which the condition U > 0 is 
an appropriate expression of the electrode effect is shown as follows. Equa- 
tion (4) is first cast into the form 


where 


Vy = ^(y,P x ) - u 


iKy,p x ) 


2Ze / A Px 
m l m 


*) 


z 2 e 2 A 2 


m 


( 6 ) 
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Therefore , 


, , 2ZeA / cp , ZeA\ 

*<y< v x> B — l Y x - i + a r) 


lim t(y,P x ) 


y-xx 


— 00 


If the function ^(y,p x ) has a form as in figure 2(a), the particle will oscil- 
late between the two points A and B as determined by the intersection of 
i)r(y,p x ) with the level of U. It is clear that a particle with U = 0 will 

reach the electrode, and, if U < 0, the particle will be absorbed or neutral- 

ized. On the other hand, if the electric and magnetic potentials have appropri- 
ate functional dependences on y, there may be such a value of p x that a nega- 
tive minimum is attained, such as at one of the points M (fig. 2(b)); then the 

particle may oscillate between the points C and D with a negative value of 
U and still not touch the electrode. (Smaller oscillation ranges are possible, 
of course, with larger values of U. ) Such cases are not suitable for use with 
the condition U > 0. When B ^ 0 is assumed, the minimum of iKy,p x ) is at- 
tained for the value of p x such that 

p x - m if = ZeA 

The extremum of is then 



and for U > 0 to he suitable, 4r m must be greater than zero. 

The distribution functions are now reformulated as functions of velocity 
rather than as functions of the integrals of motion. The following dimensionless 
forms of the variables are used: 


u^ 


v x - a 
a 





a 


u 


z 


- P 


a 


v 2 (y,v x ) 


U. 


X 


0)) 


2ZeA 

ma 


co = 


cp - clA - 
Aa 


ZeA 2 

2m 
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where 


a 2 = — 


m 


The condition U > 0 for admissible particles is also incorporated in f , and in 
terms of velocities it may be expressed as 


2ZeA 

am 


K 


cd ) - u^ r = v 2 - u 2 > 0 


or, alternatively (A < 0, Z e = -l), 


u 


x,e 


> Ofe 




,i < ^i 



(7) 


Then 



Appending the conditions on u x ensures real limits ±v for Uy. In the plasma 
at y -> oo ^ A -* -BqY, v 2 (v x ^y) -> oo , and H = 1. Since u 2 + u 2 + u 2 is bounded, 
there is the requirement that cp - <x e A and cp - a-pA both be bounded in order 
for f to be nonzero and bounded. From this requirement and from the fact that 
the electromagnetic field is independent of the species under consideration, 

a e = a i ” a 


and 


lim cp = aA + const. 


y-*co 


With the definitions 


G = G v + A (9 - OA) 
e 
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and G v defined so that 


lim G = 0 

y-H*o 


f e = C e H K,e “ “e) H ( v e 

f i = c i H K - u x,i) H ( v f - 


- u| je )exp(G)exp[-(u^ e + u^ e + u^ e )J 
ul^JexpC-ZitGjexp^ul^ + u^ ± + u^JJ 


( 8 ) 


and 



It will be shown subsequently that a is the plasma drift velocity at y = <». 
Therefore, 0 e G/e is the electric potential relative to the drifting plasma; it 
is the potential disturbance created by the presence of the electrode and is 
shown in figure 1 as the difference between the applied and the disturbed poten- 
tial curves. 


Moments of Distribution Function 

The constants C, a, and (3 may be related to plasma properties external to 
the sheath by evaluation of the moments of the velocity distribution functions. 

In the limit y 00 , H = 1, G = 0, and the plasma values of n e -> N and 
n-t -* N/Z 4 yield values for C. Also, lim uh. = 0 is obtained, from which 

- 1 ' ± y-X ) o ^ 

lim v _ = lim v Y _• = a. Then, finally, the constants 9 and 0. are related 

y-» oo x > e y-Hco x ; 1 e 1 

to the temperatures by the Boltzmann constant k by the relation 6 = kT. 

The particle densities are calculated from zero-order moments of the distri- 
bution functions as 



( 9 ) 


The other macroscopic averages are obtained, from the higher moments. Of these 
averages, the following simple results are noted.: 
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= U„ = 0 


°y = u z 


u i = \ 


From the moments the macroscopic properties are obtained: 
Mass density: 


Pe " m e n e> Pi ~ m i n i> Pm “ Pe 


+ o 




Charge density: 


a e = -en e ; a^_ = Z i en i ; a = a p + o- 


e 


Currents : 


where 


-en e v e = J e = J e + ac e i x ; J ± + aa ± l x = Z i en 1 v 1 = 


J e -ea e n e u x ^ e i x 


J;L ^i ea i n i u x,it 

j = j e + Ji = J + o-ai 


X 


J - J e + J i 

Pressure tensors for the ions and electrons are defined as 


and yield components 


p ss nma‘ 


p = nm( v - v) (v - v) 


fe - (u x) 2 ] = 2n0 [ u l - <%>’] 


p yy ‘ 2n0u y 


Pzz ™ 

Pxy = Pyx = Pzx = 0 

The off-diagonal components are zero because the velocity distribution is an even 
function of u^ and u z and the integration is symmetrical. Also, TT and u z 
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are zero for the same reason. 


SOLUTION OF EQUATIONS 
Change of Variables 

Aside from the effect of the electrode in changing the distribution func- 
tions f from the Maxwell type, the distribution functions and the densities n e 
and n-j: depend on the disturbance potential G, which is therefore a more con- 

venient variable for the problem than cp. In addition, the cutoff parameters a) 
and v depend on A in the combination given by 



where L e (= a e m e /eB 0 ) is the Larmor radius of the electrons. Although “ is 
actually a dimensionless magnetic flux, the interpretation of £ as the approxi- 
mate distance from the electrode in units of electron Larmor radius will be em- 
ployed as a more simple and descriptive concept. The system of equations (l) 
and (2) may also be reduced from fourth to third order by changing the independ- 
ent variable from y to the spatial dependence of the solution y(£) may 
then be found by a subsequent integration. The equations are 



and 


d 



B 0 L e 

BD 


where 


D , J. 5 

T 4nNe 2 c 2 

is the Debye distance, and L 2 /D 2 = 8nNm e c 2 /Bo- The method of integration of 
equation (10) is described in appendix C. 
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Partial Integration 


Equation (ll) may be integrated and the system of equations reduced to sec- 
ond order. For this purpose the momentum equation for the mixture of ions and 
electrons in a steady-state condition is utilized: 

— > — > — ► ^ 

V • (p m vv) = -V . p + oE + j xB 


vhere 


For the case at hand, 


Pm v = P 


y 

e v e 


+ P i v i 


. 

y =s i y 

e x e 


v_ 


1 


V. 
X 1 


E = iyE 
B = i z B 


a a 
” ^ = 


The relations 


V • E = 4irc^a 
VxB = 4 kj" 

inserted into the momentum equation, when integrated, yield 



vhere 

p 0 s p 0 ,i + Po,e = Ne e (i + yt) 

and Pyy is calculated from the sum of the second moments of the distribution 
functions of the ions and the electrons, that is. 
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p yy 




In order to put this result in terms of the variables being used, 

E = I a - 


/ dfi\ 

V a ' 2 dj 1 


9 9 9 

BoLe a e 
2W e = | 

SkD^c^ 


are employed, and the result is 


L /a 1 

a e dG\ 2 

9 T 2 a 2 

.. or Bg a e 

1/, . 1 \ 

n i ~ 

n e 2 

[ - (= -5 

c d Z)_ 

c 2 D 2 c 2 

2 \ + 'Ll) 

■ITV 

“TV 


(13) 


which replaces equation (ll). Because E/B is approximately equal to the drift 
velocity v x , equation (12) shows that the electric tension E 2 /8jtc 2 is always 
small compared with the magnetic pressure B 2 /8it. The magnetic and kinetic pres- 
sures are therefore always in balance, and the maximum change in magnetic pres- 
sure is approximately equal to the kinetic pressure of the plasma Pq. 


RESULTS 

The foregoing equations were integrated with consistent values for density 
and current; the results of the calculations are given in figure 3 as curves of 
the potential disturbance function G, the density ratios n e /N and nj_/N, and 
the ion and electron pressure ratios P xx /Po an( ^ ^yy/^O plotted as functions 
of £ (ssy/Lg). The pressure ratio P zz /Po is indicated to be equal to n/N. 

The calculations were made for G w = ±1 and ±5, D/L e = 100, 2, 0.5, and 0, and 
several values of cc/c and a/c. Since the effects of a, and a are very 
small, these values are not shown. In the sequence the curves for D/L e = 0, 

G w = ±1 are omitted because they are very much like the curves for D/L e =0.5. 
The method of integration used was unsuitable for the cases D/L e =0.5, G w = ±5, 
and no results were obtained. Presumably these curves would be similar to the 
curves for D/L e =0, G w = ±5. Because the coefficient of the highest derivative 
is small in these cases, standard integration techniques are not applicable. 

Also omitted is the curve for G w = -5, D/L e = 0. This carve violated condition 
(7) over a considerable region, so that the effect of the electrode is not accu- 
rately stated in condition (5). The same violation occurred for G w = -5, 

D/L e =0.5 but extended over a small region (£ < 4) where the density is low 
and presumably the effect is not large. 

Although a has an insignificant effect on the disturbance potential G, 
its effect on the actual electric potential cp relative to the electrode can be 
seen from the relations 
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f- cp = G - G w - 2 iL 5 
°e a e 


d /ecp\ dG 2a, 

“ d^ ‘ a e 

When a/a e is sufficiently large, the effect of the externally applied field 
overshadows that of the disturbance field. 


Isotropy and Equilibrium 

All the plots show some region in which either the electron or the ion pres- 
sure tensor Is not isotropic. This deviation from isotropy of the pressure ten- 
sor is a measure of the electrode effect, since it is caused by the selective re- 
moval of particles with prohibited values of the momenta p x and Py. normal to 
the magnetic force lines and by the absence of such selection on the parallel 
component p z . This same process causes the deviation of the velocity distribu- 
tion function from the equilibrium (Maxwell) distribution and reduces the number 
density below the equilibrium value Ne G for electrons or (N/Z)e _ZTG for ions. 
Consequently, the nonisotropy of the pressure tensor is considered a measure of 
deviation from equilibrium. In all cases equilibrium is established at a dis- 
tance of 4 to 5 Larmor radii appropriate to the species. 

Pressure isotropy (equilibrium) of the repelled species is enhanced by a 
electrode potential. This phenomenon results from the variation of the 
electrical field intensity (curvature of the G curve). A uniform field has no 
effect on ion collision with the anode, because a Lorentz transformation of the 
electromagnetic field to a coordinate system moving with the local drift velocity 
removes the electric field, and the ions orbit in circular paths of a size inde- 
pendent of the field strength. If the electric field is nonuniform, it cannot be 
eliminated by transformation; consequently, the orbital path is modified from the 
circular shape in such a way that repelled particles require higher energy to 
reach the electrode than in a uniform field. Thus, when the field strength is 
highly nonuniform, a very small proportion of the population of the repelled spe- 
cies is removed, and the distribution closely approximates the Maxwell distribu- 
tion with an isotropic pressure tensor. The trend to equilibrium of electrons is 
smaller than that of ions because of the reduced Larmor radius and the smaller 
variation of electric field in the orbit. These characteristics may be observed 
by comparing figures 3(a) to (d), where the sheath is calculated for various 
electrode potentials. Figures 3(e) to (h) show a similar change in approach to 
equilibrium. 


Scale of Region of Potential Variation 

The extent of the region in which the potential of electrical disturbance 
varies substantially is shown by Bertotti (ref. 4) to be of the order of the sum 
of the Debye distance and the Larmor radius of the attracted species, on the 
basis that the electrode influences the density distribution of the attracted 
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species to distances of the order of the Larmor radius, whereas the repelled spe- 
cies is in equilibrium. The distance scale of Poisson's equation is the Debye 
distance in that charge disturbances influence the potential as far as the Debye 
distance; beyond this range shielding by other charges takes place. All the cal- 
culations substantiate this conclusion (fig. 3). There is an additional smaller 
variation of potential well in a region much larger than this one in cases of 
low-voltage anodes with weak magnetic fields (figs. 3(f) and (k)), where the 
electrode potential is not sufficiently large to effect an equilibrium distribu- 
tion of the repelled species (ions); nonisotropy of the ion pressure tensor indi- 
cates this situation. 


Regions of Neutrality and Charge 

All calculations show the attainment of charge neutrality in regions farther 
from the electrode than several times the Debye length D (y !> *=5D) ; cathode 
sheaths in a weak magnetic field (figs. 3(g), (h), and (j)) exhibit only small 
charges even closer to the electrode. Poisson's equation (eq. (10)) is directly 
related to this situation in that the characteristic length for variation of po- 
tential is the Debye length; for y > 5D the derivative is small and a solution 
of the equation is the plasma condition of neutrality (n e = Zn^). 

The potential variation in the neutral region may be correlated with the 
Larmor radius and the electrode potential by consideration of equations (9) for 
density, where K represents the integral factor: 


(14) 


The integrals K are less than 1.0 by an amount that represents the fraction of 
the population that has been removed by the electrode. This interpretation then 
indicates that K -+ 1 at a distance of several times the Larmor radius since 
only the extremely high-energy particles in a Maxwell distribution could reach 
the electrode and be removed. This distance will be taken to be 5L. At smaller 
distances K -+ 1 also for a strongly repelled species. In general, > Kj_ be- 
cause of the difference in Larmor radii, but the case of an anode with a large 
potential is different in that Kp = 1 (except possibly for a very small region 
near the electrode), and in this case K e < Kj_ = 1 in the region y < 5L 0 . The 
neutrality condition 

K-j 

exp [ ( 1 + tZ)G] = T7- 
K e 

can be used with these interpretations to indicate when G > 0 (Kj_ > K e ) or 
G < 0 (K ± < K e ). 

In weak plasmas the neutral condition is attained only for y > 5D > 5Lj_. 
Neutrality thus implies K e = = 1 and, therefore, G = 0, which also implies 
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n e = Zn^ = N; that is, the plasma is undisturbed altogether in a neutral region 
(figs. 3(a) to (d)). 

If the anode potential is large (G v !> +5), ~ 1, and, therefore, when 

D < L e (dense plasma), then K e < 1 in the region 5D < y < 5L 0 . The neutrality 
condition then indicates that G > 0 (fig* 3(k)). For the region 5L e < y, 

Kg = 1 and G = 0 in any neutral region for both medium and dense plasmas 
(figs. 3(e) and (k)). 

For the cases where G v = ±1 or -5 the condition K-j_ < K e <! 1 and the neu- 
trality condition yield G < 0 where 5D < y < 5Lj_. In the range y > 5Lj_, 

= K e = 1; here the neutrality condition indicates the absence of any electrode 
disturbance at all (figs. 3(f) to (j)). 

The anode sheaths exhibit stronger charges than the corresponding cathode 
sheaths (compare fig. 3(b) with (c), (e) with (h), (f) with (g), and (i) with 
(j)). This effect is a consequence of the large Larmor radius of the ions, which 
extends the region of ion population depletion by the electrode. In a region of 
electron surplus near an anode, the ion depletion thus increases the sheath 
charge, whereas near a cathode the densities of the two species are made more 
nearly equal. The effect is large enough to effect a condition of approximate 
neutrality in the inner region (y < 5D) of the sheath in weaker magnetic fields 
(figs. 3(g), (h), and (j)) where the ion Larmor radius is larger than the Debye 
distance. 


SUMMARY OF RESULTS 

When collisions are ignored, calculations of the sheath between a plane 
electrode and a hydrogen plasma in thermal equilibrium in the presence of a mag- 
netic field parallel to the electrode show that 

1. The dimensionless disturbance potential relative to the moving plasma, 
the densities, and the pressure tensor components (ratio to undisturbed-plasma 
value) are all independent of the plasma temperature and drift velocity; they de- 
pend only on dimensionless flux (approximately the ratio of distance to Larmor 
radius), ratio of Larmor radius to Debye distance, and ratio of electrode poten- 
tial energy to thermal energy. 

2. The electrode effect on velocity distribution and pressure tensor compo- 
nents extends to a distance of five times the Larmor radius, except that the 
variation of electric field intensity (which is large with a substantially 
charged sheath) results in a reduced electrode effect on the repelled species; 
this electrode effect is less pronounced with ions near anodes than with elec- 
trons near cathodes. 

3. The plasma condition of neutrality is not satisfied in the transition 
region for rare plasmas (D > L^) nor in the region y < 5D near anodes in 
medium-density plasmas (L e < D < L^), where D is the Debye distance, L e is 
the electron Larmor radius, Lj_ is the ion Larmor radius, and y is the distance 
normal to the electrode. It is satisfied in all cases for y > 5D and approxi- 
mately so for all locations near cathodes in medium or dense plasmas (D < Lj_). 
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4. In a region of neutrality the potential is always negative in a region 
extending to approximately 51^, except that, if the anode potential is large 
(five times the thermal energy), the potential is positive in the region 
y < 5L e . 

Levis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, July 5, 1963 
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APPENDIX A 


SYMBOLS 


A 

a 

B 

b 

c 

D 

E 

e 

e 

f 

G 

g 

H 

i-x^y^z 

J 

J 

K 

k 

L 

m 

N 

n 


magnetic potential (B = -dA/dy; A(O) - 0) 

VW m - v 20 / m 

magnetic field intensity 

constant, approximation for -^/A (see eq. (C3)) 
speed of light 

Debye distance; D ~ \/kT e /4jtNe 2 c 2 
electric field intensity 
energy of particle 
proton charge 

Boltzmann distribution function 

dimensionless potential relative to plasma; G = G v + — (cp - aA) 
correction to the approximation P ~ G; g = G - P 

electrode effect function; H = 1 for permissible particles and H = 0 
for disallowed particles 

unit vectors in direction of x ; y ; or z variation 
electric current relative to plasma; J = j - ao\bc 
electric current relative to electrode 
electrode effect on density depletion (eq. (14)) 

Boltzmann constant 

Larmor radius of electrons or ions; L = a/co c 
mass of particle 

number density of electrons in plasma 
number density of particles in sheath 
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(P-H 


+ 


; )/3 


Px^y^z 


^yoz’^yy’^zz 

S 

T 

t 

U 

Sc 'V U z 


v e^i 


m 

x,y,z 

z 

a 

3 

r 

5 

€ 

5 

0 

A 

A 


xx ^yy ^zz ' 
components of generalized momentum vector 

components of pressure tensor 

right side of eq. (10) with a/e N omitted (see appendix C) 

temperature of electrons or ions 

time 

[(px + Pz)A 2 J - 2 s/ m 


nondimens ional velocity components relative to plasma; 
Ojj = (v x - a) /a; Uy = v y /a; u z = (v z - p)/a 

velocity of particle 

average velocity of electrons or ions 


velocity of ion-electron mixture; v m = 


Pe _ Pi _ 

V~ + V-t 

Pm Pm 


p e^ p i 


coordinates 

number of fundamental charge units per ion (Z 0 = -l) 
constant (drift velocity of plasma, Eq/Bq) 

constant (drift velocity in z-direction, assumed to be zero) 
approximation for G in iterative process for finding G 
-[(o/eW) + S]BqD 2 /B 2 Lg (see eq. (C3)) 

error in eq. (10) when G is replaced by T (see eq. (Cl)) 
- A / B oL e = -Ae/m e a e ; £ « y/L g 
constant (shown to be kT e or kT^) 
collision parameter (appendix B) 

-(B§L|/B 2 D 2 )[5(cT/eW)/SG + Ss/SG] (eq. (C2)) 
electrode cutoff for u y integration; v = yj^/e 2 
collision frequency per particle 
mass density of electrons or ions 
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p mass density of electron-ion mixture; p m = m e n e + m i n i 

a ’ a e’°i electric charge density of electron-ion mixture or of either component 

T T e / Tl = C J0 i 

cp electric potential function [cp( 0 ) = cp^. = 0; E = -dcp/dy] 

, 2Ze ( A Px 

* — - 

\|r stationary minimum value of 

03 electrode cutoff for u x integration; 03 = [cp - aA - ( ZeA^/2m) ] /Aa 

cjd c cyclotron frequency; 03 c = eB^/m e or 

Subscripts : 
e electrons 

i ions 

v value at electrode (y = 0) 

x,y,z component appropriate to x-, y-, or z-axis direction 
0 value in plasma (y = °°) 

Superscripts : 

( — ) velocity average 

("*) vector 


\ Z 2 e 2 A 2 

) " m 2 
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APPENDIX B 


CONDITIONS FOR WHICH COLLISIONS CAN BE IGNORED 


A limitation of the theory arises from the assumption that over a short pe- 
riod of time the effect of collisions on the distribution functions is small. 
Quantitatively this assumption is expressed by V C A> C « 1, where v is the 
scattering collision frequency and o> c is the cyclotron or Larmor frequency. 

For consideration of order of magnitude, electron shielding alone is assumed, and 
conditions are calculated in the plasma rather than in the sheath. Because of 
the large Larmor radius, ion-ion scattering gives the largest value of V C A> C - 
Then, from reference 8 v . is obtained, from which the result for hydrogen is 


where 


1 » 



A = 


5(kT) 5 / 2 
e^c^( 4:nNe^c^) 


Tj2 


The condition for ignoring collisions is 




Br 



« 0.04 



which is equivalent to a minimum magnetic field intensity. This limitation is 
very weak for high temperature and low density and very severe for low tempera- 
ture and high density. At T - 1000° K and N = 10*^ per cubic centimeter, for 
example , 
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APPENDIX C 


METHOD OF INTEGRATION 

The system of equations (10) and (ll) embodied difficulties peculiar to the 
nonlinearity of the problem, the two-point boundary conditions, and the infinite 
range of the independent variable. The method used was one of successive ap- 
proximations with linearized approximations for the corrections. With equation 
(10) written in shorter notation, 

D 2 d 2 G = _£_,o=_L + / 

" t 2 2 _„2 eN eN l 

L e B 0 d£ \ 

a function r(£) is assumed to be an approximation for G(£), and the dependence 
of S on dG/d£ is neglected “because a^/c^ and aa/c^ are small. The cor- 
rection g = G - T is assumed small , so that variations in G will cause 
changes in a and J x , which can he approximated by 

o(G,0 = a(r,c) + g(|f) G=r 


d£ a\ / a e *^x + a a \ 

2c d£ cj \ c eNa g c eNy 


j(g,0 = J(r,0 + g(^) Q=r 

If the error in the differential equation is e when r is used, 


e 




where 



(Cl) 


(C2) 


When T(0) = G w and r(») = 0, 


g (0) = g(°o) = 0 


Explicit expressions for the derivatives are 
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Because of the complicated form of the coefficient A, equation (Cl) is not 
amenable to direct solution for a self-consistent electrical field. A method of 
successive approximations was used, in which for each approximation a constant 
value of A was assumed. Then the correction g is obtained from equation (Cl) 
with the following result : 



where 

b = -/K 

and 



For L e « D and the initial approximation P = G v e convergence was very 
rapid (three iterations ) , but convergence was slow for weak fields (L e = 0(D)). 

In these cases convergence was more rapid for large values of the rate of 
convergence near the wall was improved by occasional use of values of A appro- 
priate for smaller £, although the usual procedure was to use A = (1 + Zx)L e /D, 
which is appropriate to large £. After each iteration for G, the function 
B 2 /B2 was reevaluated from equation (ll), although equation (12) might just as 
well have been used. This method was found to be unsuitable for values of 
G v > 10 or G v < -10. 

In the limiting case in which the plasma density is very large or the mag- 
netic field weak, the ratio D/L e becomes vanishingly small. Then the differen- 
tial equations (10) and (ll) are simplified by ignoring the second derivative 
terms. For this system there is the solution a = 0 or n e = Zn^, which is 
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known as the plasma equation. The plasma equation was checked by the solutions 
for the cases D/L e = l/2, G v = 1, -1 by means of integration of the differen- 
tial equation. The results showed a negligibly small net space charge. 


23 



REFERENCES 


1. Tonks, L. , and Langmuir., I.: A General Theory of the Plasma of an Arc. Phys. 

Rev., vol. 34, no. 6, Sept. 1929, pp. 876-922. 

2. Hall, Laurence S. : Probes and Magnetic Pumping in a Plasma. UCRL 6535, Univ. 

Calif. , July 19, 1961. 

3. Bobm, D. , Burhop, E. H. S., and Massey, H. S. W. : The Characteristics of 

Electrical Discharges in Magnetic Fields. A. Guthrie and R. K. Wakerling, 
eds., McGraw-Hill Book Co., Inc., 1949, ch. III. 

4. Bertotti, Bruno: Theory of an Electrostatic Probe in a Strong Magnetic Field. 

Phys. of Fluids, vol. 4, no. 8, Aug. 1961, pp. 1047-1052. 

5. Bertotti, Bruno: Theory of an Electrostatic Probe in a Strong Magnetic Field. 

II. Phys. of Fluids, vol. 5, no. 8, Aug. 1962, pp. 1010-1014. 

6. Bernstein, Ira B. , and Rabinowitz, Irving N. : Theory of Electrostatic Probes 

in a Low-Density Plasma. Phys. of Fluids, vol. 2, no. 2, Mar. -Apr. 1959, 

pp. 112-121. 

7. Landau, L. D. , and Lifshitz, E. M. : Statistical Physics. Add is on -Wes ley Pub. 

Co., 1958, p. 11. 

8. Spitzer, Lyman, Jr.: Physics of Fully Ionized Gases. Intersci. Pub., Inc., 

1956. 


24 





(Pzz/Po) 
and n~/ 



Qd/d jo n/u 


Figure 3. - Sheath structure for hydrogen plasma. Temperature of electrons same as temperature of ions 
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Figure 3. - Continued. Sheath structure for hydrogen plasma. Temperature of electrons same as temper 
ature of ions. 








Figure 3. - Continued. Sheath structure for hydrogen plasma. Temperature of electrons same as temperature of 
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Figure 3. - Continued. Sheath structure for hydrogen plasma. Temperature of electrons same as temperature of 
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Figure 3. - Continued. Sheath structure for hydrogen plasma. Temperature of electrons 
temperature of ions. 
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Figure 3. - Continued. Sheath structure for hydrogen plasma. Temperature of electrons same as temperature 
of ions. 








Figure 3- - Continued. Sheath structure for hydrogen plasma. Temperature of electrons same as temperature of 














Figure 3. - Continued. Sheath structure for hydrogen plasma. Temperature of electrons same as temperature of 




Figure 5. - Continued. Sheath structure for hydrogen plasma. Temperature of electrons same as temperature of 










